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COPAR - Multivariate time series modeling using 
the COPula AutoRegressive model 
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Abstract 



Analysis of multivariate time series is a common problem in areas like finance and eco- 
nomics. The classical tool for this purpose are vector autoregressive models. These how- 
ever are limited to the modeling of linear and symmetric dependence. We propose a novel 
copula-based model which allows for non-linear and asymmetric modeling of serial as well 
as between-series dependencies. The model exploits the flexibility of vine copulas which are 
built up by bivariate copulas only. We describe statistical inference techniques for the new 
model and demonstrate its usefulness in three relevant applications: We analyze time series 
of macroeconomic indicators, of electricity load demands and of bond portfolio returns. 
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00 1 Introduction 

(N 

£<■} The analysis of multiple time series is of fundamental interest in finance and economics. Clas- 

^X sically, interdependencies among multivariate time series have been modeled using vector au- 

toregressive (VAR) models. Such models provide insights into the dynamic relationship of the 
time series and often produce forecasts superior to independent univariate models. VAR models 
in economics were advocated by |Sims (1980)[ standard reference books are Liitkepohl (2005) 



Hamilton (1994)| and [Tsay (2002) 



The bivariate pth order vector autoregressive model, VAR(p), for two time series {Xt} and 



{Y t } is defined as 





+*! :r +... +*,{?-*) + {:) . (i.i) 



where 3>j, j = 1, ...,p, are 2-by-2-matrices of autoregressive coefficients and c\ and C2 are con- 
stants. The vector et = (£±,£2)' is multivariate white noise, that is E{et) = and E(etS s ) = E 
for t = s and otherwise, where E is a symmetric positive definite 2-by-2-matrix. Typically 
Et ~ ^2(0, E) is assumed. 
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While VAR models can only capture linear and symmetric dependence in time and between 
series, we propose a new copula-based model which overcomes such limitations and allows for an 
extremely flexible modeling. Copulas are the canonical statistical tool for statistical dependence 



modeling. The theorem by Sklar (1959) shows that every multivariate distribution can be 
represented in terms of a copula which couples the univariate marginal distributions. For a 
random vector X = (X\, ..., X d ) ~ F with marginal distributions Fj, i = 1, ..., d, it is 

F( Xl , ..., x d ) = C(F 1 (x 1 ), .., F d (x d )), (1.2) 

where C is some appropriate (i-dimensional copula, a multivariate distribution on the unit hy- 
percube with uniform margins (see Nelsen (2006)|and Joe (1997) for more details). 



We will use a fully integrated copula model to capture effects in time and between series. In 



particular, our model is built upon a so-called vine copula (see Kurowicka and Joe (2011) for an 



overview) . Such vine copulas are flexible multivariate copulas constructed through a sequence of 



bivariate copulas, a pair-copula decomposition. While Smith et al. (2010) recently showed how 
univariate time series can be modeled using a so-called D-vine pair-copula decomposition, we will 
show how such pair-copula decompositions can be conveniently used to model the dependence 
among multiple time series. 

The contributions of this paper are as follows: We introduce the so-called copula autoregres- 
sive model, COPAR, which exploits the enormous flexibility of vine copula models and allows 
for non-linear and non-symmetric modeling of serial and between-series dependence. By allow- 
ing for arbitrary marginal distributions, the model can also account for common features of 
univariate economic and financial time series like skewness and heavy-tailedness which are not 
captured appropriately using a normal distribution. Required statistical inference techniques for 
the model are presented and described in detail. In addition, we also discuss how the model can 
be easily used to test for Granger causality, a central concept to determine interdependencies 
among multiple time series. The usefulness of our model is demonstrated and carefully evaluated 
in three relevant applications: We analyze monthly macro-economic indicators, daily electricity 
load demands as well as monthly bond portfolio returns. 

The paper is structured as follows. In Section [2] we establish the relevant technical back- 
ground on copulas and pair-copula constructions in particular. The copula autoregressive model 
is introduced and discussed in detail in Section [3j The three applications are subsequently 
treated in Section |4j while Section [5] concludes with an outlook to future research. 

2 Pair-copula decompositions for univariate time series 

Let {Xt}t=i....,T be a univariate time series of continuously distributed data. Through condi- 
tioning the joint distribution of {Xt} can be decomposed as 

T 

f(x X ,...,X T ) = f{xi)Y[ft\l : (t-l)(xt\xi,-,Xt-l), (2.1) 

t=2 



where / denotes the common marginal density of Xt, t = 1, ...,T; F will denote the correspond- 
ing distribution function. Here we use r : s := (r, r + 1, ..., s — 1, s) for r < s and / s m denotes 
the conditional density of X s given {X r ,r £ D}. 



As outlined in Smith et al. (2010) this expression can be used to obtain a general de- 
composition in terms of bivariate copulas. For the distribution of X s and Xt, s < t, given 
{X s+ x, ...,X t -i} (to shorten notation we often write X( s+ xyui\) it follows according to Sklar's 



theorem (1.2) that 



fs,t\(s+l):(t-l)(x s ,X t \x s+ l,...X t ^i) = 

Cs,t\(s+l):(t-l){F S \(s+l):(t-l)( x s\x s +l, ■■■, X t -l), F t ^ s+1 y(t-l)(xt\x s +l, ..., Xt-l)) 
X fs\(s+l):(t-l)(Xs\x s +l,---,X t -l) X /t[( a +l):(t-l) (x t \x s +l, ■■; Xt-l)) 

where c g t|( s +i) : (t-i) is an appropriate bivariate copula density. Rearranging terms gives 

ft\a:(t-l)(xt\Xs, -Xt-l) = 

C S ,t\(s+l):(t-l)(Fs\(s+l):(t~l)(Xs\x s +l, ...,X t -l), F t \( s+1 y.(t-l)(x t \x s+ l, —,%t-l)) 
X ft\( S +l):(t-l)(xt\x s+ i,...,X t -i). 



Let u r |( r+ i).(t_i) := i ? r .|( r . +1 ).( t _ 1 )(x r |x. r+ i, ..., 

t- 



xt-i) andn t | (r . + i).( t _i) := ^(^.(^(xtlav+i, ...,x t -i). 
1 one then obtains 



By recursive conditioning on s = 1,2, 
ft\i-.(t-l)(xt\%i, ...x t -i) = 

f{xt)ct-i,t{F{x t -i),F(xt)) \\ c rjt |( r . +1 ).( t _ 1 )(u r |( r+1) .( t _ 1 ),« t |( r . +1 ).( t _ 1 )) 



t-2 



(2.2) 



r=l 



Finally, plugging Equation (2.2) into Expression (2.1) gives 

T T t-1 

/(xi,...,Xt) = 11/(^)11 \_\_ c r,t\(r+l):(t-l)( u r\(r+l):(t-l)i u t\(r+iy.(t-l))- 
t=\ t=2r=l 

This is the product of T(T — l)/2 bivariate copulas, so-called pair-copulas, and the marginal 
densities evaluated at each time point Xt, t = 1,...,T. The construction does not require the 
selection of a particular copula family, so that very flexible models can be deduced from it. 
However, it is clear that copulas corresponding to the same time lag have to be identical. For 
example Cj_2,t-i and Ct-\,t must not be different. 

The above construction is called pair-copula decomposition and was introduced by|Aas et al 



(2009) The particular way described here is called D-vine and belongs to the more general 



class of regular vines (R-vines) introduced by Joe (1996) and Bedford and Cooke (2001, 2002) 



and described in more detail in jKurowicka and Cooke (2006)| and Kurowicka and Joe (2011 



R-vines are a graphic theoretic model to determine which pairs are included in a pair-copula 



decomposition. The following definition is taken from Kurowicka and Cooke (2006) 



Definition 2.1 (Regular vine). A regular vine (R-vine) on d variables is a sequence of linked 
trees (connected acyclic graphs) 71, ..., Td-\ with nodes M% and edges Si for i = 1, ..., d— 1 which 
satisfy the following three conditions. 



Xx,X 2 x—v x 2 ,x 3 ^-v X 3 ,X± /-—-. x 4 ,x 5 
Xi ) — (X 2 ) — (X 3 ) — (X 4 ) — [X 6 



Ai , A3 J X2 ^**~ " -»«^ *\.% 1 A4 J .A3 ^*~~ " — ■ .^^ A3 , A5 J A4 

X^Xi ) Q X 2 ,X 3 ) C -^3,-^4 ) C x^,x 5 



c2^ 



-^4|^2:3 J C ^2^5^3:4 



Figure 1: D-vine for T = 5 with edge labels. 



(i) Tree 7i has nodes A/i = { 1 , . . . , d}. 
(ii) For i = 2, ..., d — 1 tree 7i has nodes A/i = £j_i. 

(iii) If two edges in tree % are joined in tree 7i+i, they must share a common node in tree %. 
The last property is called proximity condition. 

D-vines are R-vines, where each node is connected to at most two other nodes. The D-vine 
corresponding to the above decomposition is shown in Figure [TJ 

By associating each edge e = j(e),k(e)\D(e) in an R-vine with a bivariate copula density 
c j(e),k(e)\D(e)i the complete pair-copula decomposition is defined. The nodes j(e) and k(e) are 
called conditioned nodes and D(e) the conditioning set, where in each tree from top to bottom 
an additional variable is added in the conditioning set of the bivariate copula. 

Theorem 2.2 (R-vine density (Kurowicka and Cooke 2006, Theorem 4.2)). The joint density 
of X\ , ... , Xd is uniquely determined and given by 



f(xi,...,x d ) 



Y[fk(xk) 



.k=l 



d-1 



II II C j(e),k(e)\D(e)(. F ( x j(e)\ x D(e))i F ( x k(e)\ x D(e))) 



(2.3) 



where x D ^ denotes the sub-vector of x = (xi, ...,Xd)' determined by the indices in D(e). 



In (2.3) the arguments of copulas in tree T% can be recursively computed from copulas in 



trees 7i, . . . , 77_ 1 using the general formula 

9C XV . 1 „ _ ( F ( x I v _ j ) , F (vj I v - j ) ) 



F(x\v) 



dF(vj\v 



(2.4) 



3\ v ~3j 



where C xv .\ v _. is a bivariate copula, Vj is an arbitrary component of v and v^.j denotes the 
vector v excluding Vj. 



To facilitate statistical inference of R-vines, they can be conveniently stored in matrix no- 



tation as recently proposed by Morales-Napoles (2011) and further explored by Difimann et al 
(2012) Let M £ {0, ...,d} be a lower triangular matrix, where the diagonal entries of M 
are the numbers 1, ..., d in decreasing order. In this matrix, according to technical conditions, 
each row from the bottom up represents a tree, where the conditioned set is identified by a 
diagonal entry and by the corresponding column entry of the row under consideration, while the 
conditioning set is given by the column entries below this row. Corresponding copula types and 
parameters can conveniently be stored in matrices related to M. The fixed ordering of diagonal 
entries ensures uniqueness of the R-vine matrix. 

The serial D-vine decomposition presented above can be stored in the following matrix 

\ 



/ X T 




Xi 


Xt-i 


x 2 


X! 


Xt-2 


Xt-3 


\ Xt-i 


Xt-2 



x 3 

X\ X2 
X2 X\ 



x 1 ) 



which is easily extendible to include future observations {Xt+i,Xt+2, ■•■}■ For example, the 

second entry in the first column identifies the conditioned pair X\ and Xt given {X2, ..., Xt-i}- 

Corresponding copula types are stored in the off-diagonal entry associated with the pair: 



/ 



\ 



C 



1,T|2:(T-1) 



C2,T|3:(T-1) W,T-l|2:(T-2) 



Ct-2,T\T-1 Ct-3,T-1|T-2 



C 



1,3|2 



V Or- 



1,T 



Ct-2,T-1 



O 



2.3 



c 



1,2 



/ 



Here in each row the same copula type must be used. For deriving the joint likelihood the 
pair-copulas have to be evaluated in conjunction with the conditional distribution functions. 

Using this matrix notation Difimann et al. (2012)| give algorithms to compute the log- 
likelihood of an R-vine and to sample from an R-vine. Thus maximum likelihood estimation 
of copula parameters is feasible. Copula types are often selected sequentially starting from the 



first tree (see also Brechmann et al. (2012)). 



3 The copula autoregressive model 

Let {Xt}t=i,...,T and {lt}t=l,...,T be two univariate time series jointly observed at time point 
t = 1, ...,T. The aim of this paper is to derive a flexible multivariate distribution of {X t } and 



{Yt}, in particular allowing for non-linear dependence — serial as well as between-series. On the 
one hand, such a distribution can be used to investigate the dependence among the time series 
(in-sample fit). This involves for instance testing for Granger causality which will be discussed 
below. On the other hand and most importantly, future values can be predicted based on this 
distribution (out-of-sample prediction). 

Our model is based on a particular R-vine structure and defined as follows. 

Definition 3.1 (COPAR). The copula autoregressive model (COPAR) for time series {Xt}t=i,...j 

and {Yt\t=\,...,T has the following components. 

(i) Unconditional marginal distributions Fx and Fy of {Xt} and {Yt}, respectively. 

(ii) An R-vine for the serial and between-series dependence of {Xt} and {Yt}, where the 
following pairs are selected. 

(a) Serial dependence of {X t }: The pairs of a serial D-vine for X\, ...,Xt, i.e. 

X s ,X t \X s+1 ,...,X t - 1 , \<s<t<T. (3.1) 

(b) Between-series dependence: 

X s ,Y t \X s+1 ,...,X t , l<s<t<T, (3.2) 

and 

y a ,x t |Xi,...,x t _i,y a+1 ,...,y t _i, i< s <t<T. (3.3) 

(c) Conditional serial dependence of {Yt}: The pairs of a serial D-vine for Y\,...,Yt 
conditioned on all previous values of {Xt}, i.e. 

y SJ yt|Xi J „. J x t ,y B+ i,... ) y 4 _i ! i< s <t<r. (3.4) 

Pair-copulas of the same lag length t — s, t > s, are identical. We associate 



(a) copula Cf_ s := C XsA | X with Expression pi] ), 



(b) copulas Cg? s := C Xs ,Y t \x {s+x) .. t and Cj_ x s := C Y „x t \x Ht _ 1 - i ,Y(. +mt -i) with Expres- 



sions (3.2) and (3.3), respectively, and 



(c) copula C(_ s := C Ys ,Y t \Xt :t ,Y (s+1) .. {t _ x) with Expression pT4). 



Remark 3.2. The number of different pair-copulas utilized in the COPAR model of Definition 
3.1 is AT — 3: For C£_ s , Cjsg and Cj_ s with s < t there are T — l different ones each. In addition, 
there are T different copulas C^X, since s < t. 



The joint density of the COPAR model can be derived through Expression ( 2.3 ) . To illustrate 
the rather technical definition, we present a small-dimensional example. 



Example 3.3. Let T = 4. Then the COPAR model of Definition 3.1 for the variables X 1 ,...,X 4 



and Y\, ..., Y± is constructed as shown in Figure [2j □ 




(3.1) 




<^x"i,x 4 |x 



CYuY 2 I X 



(3.3) 




-C^X 3 \X 



C^YuX 3 \X 1[2 ,Y 2 



(3.4) 



<^y 3 |-x: 



(3.3) 



C^Y^X^^y 



(3.3 



(3.4) 




< ^,x 4 |x^ r^ — — — (^x i,y 4 \x7 t^) 

(3-4) _ — _ 

C\,Y 4 \Xi£} 



Figure 2: Example for T = 4. Edge labels relate to Definition 3.1 



Clearly, {Xt} plays a pivotal role in this modeling approach: While the serial dependence 
of {Xt} is modeled unconditionally, that of {Yt} is specified conditionally on {^t}. In other 
words, the roles of {Xt} and {Yt} cannot be simply interchanged, since the time series play 
different roles. On the other hand, the R-vine structure specifies the full joint distribution of 
{Xt} and {Yt} so that this is mainly an issue of interpretability. Its implications will be treated 



in more detail when prediction is discussed (see Section 3.3). However note that the fitted joint 
distribution may in the end actually look differently depending on the order of variables in the 
modeling, since different copulas can be used. In the case of Gaussian pair-copulas, this is not 
the case, since then the R-vine copula corresponds to a multivariate Gaussian copula, where the 



correlation matrix can be computed from conditional correlations as given by the R-vine copula 
parameters. 

While the serial dependence is rather straightforward to understand from this model, the 
modeling of the between-series dependence warrants a more detailed examination. For this 
purpose it is useful to look at the R-vine matrices associated to the R-vine structure of Definition 



3.1 Here, we continue with Example 3.3 first. 



Example 3.4. The R-vine matrix corresponding to the R-vine structure in Figure [2] can be 
derived as 

\ 



( YA 










Yi 
Y 2 
Y 3 

x x 

x 2 

x 3 

\ x 4 


X A 




Y x 
Y 2 
Y 3 


Y 3 




Y 2 


x 3 




Y 2 


Y 2 


x x 

x 2 

X 3 


Xi 

x 2 
x 3 


Y X 


X x 

x 2 


X x 

x 2 



x 2 

Yx 



Yx 



X, X, 



(3.5) 



Xx) 



The solid line is drawn to highlight the structure. On the diagonal, values of {X{\ and {Yt} 
appear alternately starting with X\ and increasing from right to left. It is clear that in this way 
the matrix can easily be extended to include new observations as additional columns on the left. 



This will prove useful for forecasting as discussed in Section 3.3 To gain detailed insight into 
the structure and the dependence properties the model implies, we take all odd numbered and 
all even numbered columns, that is columns 1, 3, 5 and 7 and 2, 4, 6 and 8, respectively, and 
look at them separately. Note however that these sub-matrices are not valid R-vine matrices 
themselves. 

• Even numbered columns: 



X A 






Yx 






Y 2 


x 3 




Y 3 


Yx 




Xx 


Y 2 


x 2 


x 2 


Xx 


Yx 



(3.6) 



\ X 3 X 2 Xx | Xx J 
Now it becomes clear that the pairs below the solid line specify the serial dependence of 



{Xt} (see Expression ( |3.1| )). The pairs above the solid line however model between-series 
dependence: Given past values of {X t } (and of {it}), dependence of {X t } with regard 



to previous values of {Y t } is modeled (see Expression (3.3)). The first column gives the 
following pairs: F 3 ,Xi|X 1: 3, Y 2 ,X 4 \Xx- 3 ,Y 3 and Yx,X 4 \Xx- 3 ,Y 2:3 . 



Odd numbered columns: 



( n 




\ 


Y 1 






Y 2 


Y 3 




Y 3 


Yy 




x 2 


Y 2 


Y 2 


X X 


Yi 


x 3 


x 2 


Xi Yt. 



V x 4 



^3 



(3.7) 



X 2 Xi ) 



Here, the pairs below the solid line also model between-series dependence, namely that of 
{Yt} with respect to previous values of {Xt} given realizations of {Xt} between the two 
time indices under consideration (see Expression ( |3.2[ )). For example the first column spec- 
ifies the following pairs: X&,Yi (unconditional dependence), X^jY^Xi, X 2 ,Yi\Xs-A and 
Xi,Yi\X 2: /i. Finally the pairs above the solid line specify the serial dependence structure 
of {Yt} conditioned on all observed values of {Xt} up to the maximal time index of the 



pair (see Expression (3.4)). If {X t } and {Y t } are independent, the dependence structure 



of {Y^} also is an unconditional serial D-vine. 
The R-vine copula matrices corresponding to these two sub-matrices are then given as follows: 



Even numbered columns: 
( 



\ ( 



a 



Y 1 X 4 \X 1:3 Y 2 :3 



c 



Y 2 X 4 \X 1 .. 3 Y 3 



c 



y 3 x 4 |Xi;3 



a 



c 



X\X 4 \X 2 



Y X X S \Xl: 2 Y 2 



a 



Y 2 X 3 \X 1: - 



Cx 2 x 4 \x 3 
V Cx 3 x 4 



c 



X\X 3 \X 2 



a 



YjX2\Xi 



Cx 



2X3 



CxiX 2 



cl x 

YX 



c\ 



c x 
c x 



cl x 

C\ x 



c 2 x 

C x 



C YX 



• Odd numbered columns: 
( 



a 



\ 



Y 1 Y 4 \X 1 . A Y 2 :3 



a 



Y 2 Y 4 \XvaY 3 



c 



Y3Y4XU4 



c 

c 



X\Y 4 \X 2 - 4 
X 2 Y4,\X 3 -A 



c 



X 3 Y 4 \X 4 



Cx 



4 Y 4 



c 



Y 1 Y 3 \X 1 . 3 Y 2 



C 



Y 2 Y 3 \Xj 



c 



XlY 3 |X 2 :3 



c 



X 2 Y 3 \X 3 
Cx 3 Y 3 



c 



y 1 y 2 \x 1 . 



c 



x x y 2 \x 2 

Cx 2 Y 2 Cx 1 Y 1 J 



cl 

cX 



c\ 


cl 
cl 




c? Y 




cr 


c XY 

C XY 


c\ 




C XY 


C XY 





\c£ 



XY 



r>XY 



r<XY 



r>XY 



This illustrates the notation of pair-copulas in Definition 3.1 For the serial D-vine dependence 
of {Xt} and {Yt} it is clear that copulas for the same lag length are identical. For instance 



X t - t 

p- 

Yt-i 



%X t 

,y 

t 

*Y t 



Figure 3: Simplified illustration of the dependence structure of the COPAR model. Solid lines 
represent serial dependence, dashed and dotted lines between-series dependence given by the 



entries specified below the solid line in matrix (3.7) and above the solid line in matrix (3.6) 
respectively. 



Cf = Cx 1 x 2 = Cx 2 x z = Cx 3 x 4 - This translates to the other copulas: Copulas with the 
same lag length of the conditioned nodes are identical, e.g. C\ x 



a 



Y 1 X 2 \X 1 



c 



Y 2 X 3 \X 1: 2 



a 



YsXAXw When combining these two matrices, one obtains the complete R-vine copula matrix 
corresponding to the structure matrix (3.5) utilizing 13 different pair-copulas. □ 



Similar to Matrix (3.5) the general R-vine structure matrix of the joint distribution of 
{Xt}t=i,...,T and {Yt}t=i,...,T as defined in Definition 3.1 can be derived as 



/ Y T 
Yi 

Y 2 



\ 



Xj< 
Yl 

Y, 



Y T - 



i 



Xx 
X 2 



Y T -i 
X x 
X 2 



Y T -i 
Yl 
Y 2 

Y T - 2 
Xi 
X 2 



Xt-i 

Yi 
Y 2 

Y T -i 



X x 

X 2 



\ JCt Xt—i Xt—i Xt— 2 



X 3 








Y x 


Y 2 




Y 2 


Yx 


x 2 




Xx 


Xx 

x 2 


Yx 


Yx 


x 2 


Xx 


Xx | 



(3.8) 



Xx J 



The interpretation of this 2T-by-2T-matrix is the same as in the above example. The corre- 
sponding copula matrix is also found in exactly the same way and therefore not shown here. 
Figure [3] depicts the dependence structure specified by the different blocks of the matrices. 



As noted in Remark 3.2 the COPAR model requires AT — 3 copulas for dependence modeling, 
while a standard R-vine specification would use T(2T— 1). In other words, the number of copulas 
grows only linearly in the number of time points in contrast to quadratic growth of a standard 
R-vine. This number even reduces when assuming an autoregressive structure of a specific order. 
We therefore define the following model. 



10 



Definition 3.5 (COPAR(/c)). The copula autoregressive model of order k (COPAR(fc)) is de- 
fined as the COPAR model of Definition [3TTJ where all pair-copulas corresponding to a lag length 
greater than k are independence copulas, that is 



a 



X 



a, 



XY 

t-s 



a 



YX 



a 



t-s 



II for t — s > k, 



where II denotes the independence copula. 



Remark 3.6. Let k denote the autoregressive order, then the number of different pair-copulas 



in the COPAR(/c) model of Definition 3.5 is 4k + 1. 



Clearly, the COPAR model of Definition 3.1 then corresponds to a COPAR(T — 1) model. 



We continue Examples 3.3 and|3.4|to illustrate this definition. 



Example 3.7. If we assume a COPAR(l) model for the dependence of X\, ...,X% and Y\, ..., Y4, 



some of the copulas specified in the R-vine copula matrix of Example 3.3 simplify to indepen- 
dence copulas. For example C Y and C XY are independence copulas. The sub-matrices that we 
considered before can therefore be reduced as follows, where independence copulas are indicated 
by the symbol II. The left matrix corresponds to even numbered columns, the right to odd 
numbered columns: 



/ 



\ 



/ 



n 
n 



c YX 


n 


n 


c YX 


n 


n 



V c? c? 



x 



C{ 



I 



\ 



n 
n 


n 


C Y 




n 




n 


n 

c XY 

r<XY 




c xy 


C XY 
r<XY 




cr 


cf Y 



The number of different pair-copulas thus reduces from 13 to 5. Similarly, assuming an autore- 
gressive structure of order 2 would still reduce the number of copulas to 9. □ 



Clearly, the number of different pair-copulas (see Remark 3.6) no longer depends on the 
number of time points T and thus allows for a very parsimonious modeling. For example, if an 
autoregressive structure of order 2 is assumed, only 9 different copulas are required for describing 
the dependence of {Xt}t=i,. t an d {Yt}t=i,...,T, while T may be 1000 or greater. In terms of 
model parameters this means that only a small number is needed, since most common bivariate 
copulas have at most two parameters and also marginal distributions have rarely more than 
three or four parameters. In particular, compared to the model complexity of the bivariate VAR 



model (1.1) there is essentially no difference: While a VAR(fc) model requires 4k + 1 parameters 
(k 2-by-2-matrices with 4 entries each as well as the off-diagonal entry S12 of the covariance 
matrix E) for modeling the serial and the between-series dependence, the COPAR(/c) model 
needs exactly the same number if only one parameter copulas are used. Additionally, a VAR 
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model uses two constant regressors as well as two residual variances (the diagonal entries En 
and £22 °f the covariance matrix S) to model the margins. In contrast, a COPAR model needs 
the full specification of marginal distributions, which however usually does not require many 
parameters (e.g. for skew- normal margins as used in Section HI below, six parameters are needed 
for the COPAR model in contrast to four needed by the VAR model). 



3.1 Extension to higher dimensions 

The above construction for two time series {X t }t=x,...,T and {Yt\t=i,.,.,T can also be extended 
to an arbitrary number of time series. The following example illustrates the case of three time 
series {X t }, {Y t } and {Z t }. 



an 



Example 3.8. Again let T = 4. Similarly to Definition 3.1 and Examples 3.3 and 3.4 
appropriate R-vine structure for the variables X\, ...,Xx,Yx, ...,Yx, Zx, ..., Z4 is constructed as 
the following R-vine structure matrix: 
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(3.9) 



y Xx x A X 3 X 3 X 3 X 2 X 2 X 2 Xx Xx Xx I Xx / 

The patterns clearly resemble those of the two-dimensional case. Serial dependence of {Xt} 
and {Yt} is again modeled using serial D-vine structures (for {Yt} conditionally on observed 
values of {X t }). In the same way, the serial dependence of {Z t } is captured by a serial D-vine 
structure conditioned on observed values of {X t } and {Yt}, that means in terms of the pairs 
Z s ,Zt\Xi,...,Xt,Yi,...,Yt,Zg+i,...,Zt-i, 1 < s < t < 4. Between-series dependence of {Yt} and 
{Z t } is also specified conditionally on values of {X t }. □ 

Along the lines of this example, multivariate time series can be modeled by iteratively con- 
ditioned D-vines and appropriate between-series copulas. Let m E N be the number of differ- 
ent time series {X t j}, j = l,...,m, where {Aji} is the pivotal time series, {X t2 } the second 
pivot, and so on. Using the notation of Definition |3.1[ an adequate R-vine based autoregres- 



sive model can then be specified in terms of m blocks of T — 1 pair-copulas C t _ J s each for 



(conditional) serial dependence of {X t j}, j 



1, ...,m, as well as 



(?) 



blocks of T pair-copulas 



Ct-s' LJ each for between-series dependence of {Xu} and {Xtj} with i < j and (™) blocks of 
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T — 1 pair-copulas C t J s * each for between-series dependence when i > j. The number of pair- 
copulas used in such a model is m(T — 1) + (™)T + (™) (T — 1) = m 2 T — m(m + l)/2 and can 
again be significantly reduced by assuming an appropriate autoregressive order k, namely to 
mk + (™) (k + 1) + (™)fc = m 2 k + m(m — l)/2, which is also the number of parameters used in 
a VAR(/c) model for between-series dependence of m time series (not counting parameters used 
for marginal modeling). 

This modeling approach opens up new possibilities in constructing flexible autoregressive 
models for arbitrary numbers of time series. For simplicity and for illustrative reasons we 
concentrate here on the case of two time series. 

3.2 Model estimation and selection 

In this section we discuss several issues related to selection and estimation of COPAR models. 
We begin with a note on the marginal distributions. 

The problem about the selection of appropriate marginal distributions is that no i.i.d. obser- 
vations are available based on which the distribution could be chosen. The selection is therefore 
subject to uncertainty and should be handled carefully. We hence advocate using rather more 
complex versions of common univariate distributions, such as skew-normal or skew-t distribu- 
tions, to be able to capture features of the data which may be otherwise disguised through the 
serial dependence. 

With respect to parameter estimation, all parameters of a COPAR(A;) model can be estimated 
jointly by maximum likelihood techniques, since the model is rather parsimonious. A common 



alternative in the literature is estimation using inference functions for margins (IFM) by Joe 



and Xu (1996)| This means that first the parameters of the marginal distributions are estimated 



and then, given these parameters, the copula parameters. As noted above, there is no i.i.d. data 
available for the margins, so that IFM estimation is not possible in our case. Given a good 
selection of the marginal distributions, IFM-type estimation ignoring the serial dependence is 
however useful for selection of copula types in a pre-parameter-estimation step. 

When constructing a COPAR(/c) model, Ak + 1 different copulas need to be chosen, where k 
denotes the autoregressive order. Since for example the copula Cx 1 x 3 \x 2 depends on the copulas 
Cx±x 2 an d Cx 2 x 3 through its input arguments, copula selection will proceed sequentially. Figure 
H^ depicts the interdependencies of the copulas for a COPAR(2) model, more details are given in 
Appendix [Aj Copula selection itself can be done for example using the AIC to penalize copula 
families with more parameters. Note that models can also be estimated in this sequential way. 
Resulting IFM-type estimates are typically good starting values for a full maximum likelihood 
estimation. 

Finally, the autoregressive order has to be selected. We propose two methods. First, when 
selecting copulas, an independence test can be used to check whether the independence copula 
is appropriate. If all copulas with lag greater than k* are selected as independence, then k* 
is the selected autoregressive order of the COPAR model. Alternatively, one may fit different 
COPAR(/c) models for different lag lengths k > 1. The optimal lag k* can then be chosen 
such that the COPAR(A;*) model minimizes an information criterion such as the Akaike (AIC), 
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Figure 4: Interdependencies of pair-copulas in a C0PAR(2) model. 



Bayesian (BIC) or Hannan-Quinn (HQC): 

AIC(/t) = - 

BIC(/t) = - 

HQC(fc) = - 



24 + 2pk, 

2£ k + log(2T) Pk , 

24 + 21og(log(2T))p fe , 



(3.10) 



where 4 denotes the estimated log likelihood of the COPAR(/c) model and p k the number of its 
parameters (4/c + 1 plus marginal parameters). 



3.3 Forecasting 

A major purpose of autoregressive modeling is forecasting. The autoregressive R-vine model 
can easily be used for this. Given time series {X t }t=i,...,T and {Yt\t=l,...,T, we like to forecast 
XT+h and Yr+h, where h > 1 (/i-step prediction). 

In the case h = 1 (one-step prediction) this can be established iteratively using the following 
decomposition of the distribution function of Xt+i, Yr+i\Xi, ...,Xt, Yy, ...,Yp: 

F{x T+ i,y T+1 \xi, ...,x T ,yi, ...,y T ) = F{x T+1 \xi, ...,x T , yx, ...,y T )F(y T+1 \xi, ...,x T +i,yi, -,yr), 

where the two univariate conditional distribution functions can be described in closed form using 
pair-copulas of our R-vine model. By conditional inverse sampling first of Xt+i\Xi, ...,Xt, Y±, ..., Yt 
and then of Yt+i\Xi, ..., Xt+i, Yi, ■■■,Yt, a forecast can be derived. This means that the selec- 
tion which time series corresponds to {Xt} and {It}, respectively, determines which variable 
can be directly predicted and which conditionally. 

If h > 1, then X^+h and Yr+h can be predicted in essentially the same way by iteration: 
first predict Xt+i and Yp+i, then Xt+2 and Yp+2 and so on. 

An illustrative example provides more details. 



Example 3.9. In the setting of Example 3.3 we would like to predict X§ and I5 given X-y, ..., X4 



and Yi, ..., I4. Since the dependence model of the latter eight variables is already known, the 
additional variables X5 and I5 have to be integrated into this model appropriately such that 
we are able to determine the conditional distribution of X^jY^lXx, ...,X^ Y\, ...,Yi. This is 



straightforward using the model building principles of Definition 3.1 and illustrated in Figure 
10 in Appendix [Bj 
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In terms of the R-vine structure matrix this means the addition of two new columns: 
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Using Expression 2.4 the conditional distribution of Xs\Xx, ..., X4, Yi, ..., I4 can then be derived 



from this structure matrix as 

dC Yl x 5 \x 1 . A Y 2 . A ( F ( x 5\ x i, ■■■, x±, V2, •••, 2/4), F(yx\xx, ..., x 4 , y 2 , ..., y 4 )) 



F(x 5 \xx,...,Xi,yx,...,yi) 



dF(yx\xx,...,Xi,y 2 ,...,yi) 



where the arguments can be iteratively decomposed in the same way in terms of bivariate copulas 
specified in the model. For instance, 



F(yx\xx,..-,Xi,y 2 ,...,yi) 



dCvxY^XuiYwiFiyAxx, ..., x^ y 2 ,y 3 ),F(yi\xx, ..., x A , y 2 , y 3 )) 



dF(yi\xx,.-,Xi,y 2 ,y 3 ) 



and so on. However note that the model requires copulas C, 



x 



CxiX 5 \X V a and Q 



Cy 1 x 5 \x 1A Y 2A - These copulas however are not known from the model and have to be taken as 
independence copulas. In general, this will never be a problem, since in a COPAR(/c) model 
with k < T — 1 (here T = 4) these pair-copulas are independence copulas anyway and a forecast 
based on only four values is not quite sensible but only used for illustration here. 

The case of the distribution Y$\Xx, ■■■, X^jYx, ...,Yi is similar and therefore not discussed 
here in more detail. □ 

Clearly, forecasting Xx+h and Yp+h given Xx, ...,Xy, Yx, ■■■,Yx simply requires to extend 
the R-vine matrices according to the model construction principles and then determining the 
conditional distribution functions in terms of bivariate conditional distribution functions which 
are derived from the corresponding bivariate copulas. An appropriate autoregressive structure 
of order k ensures that all required copulas are well-defined and also simplifies the computational 
burden significantly, since not all observed values have to be taken into account in the modeling. 

Finally a note on conditional forecasting. In order to forecast Yr+i given Xx, ...,Xt, Yx, ...,Yt, 
we derived the conditional distribution of Yt+x\Xx, ■■■,Xt+x,Yx, ...,Yt, which depends on the 
value of Xt+x- To obtain a joint forecast of Xt+x and 1t+i, the variable Xt+x is also predicted. 
If however the true (observed) value of Xt+x is employed here, this is called conditional forecast- 
ing, which is a very attractive tool of autoregressive models. On the one hand, it is interesting 
in economics, since, for example, economic indicators such as inflation and unemployment rates 
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are often not released simultaneously. Let us assume that inflation rates are released first. This 
information could then be used to obtain more accurate forecasts of unemployment rates. Ad- 
ditionally, conditional forecasting is very useful for scenario analysis, for instance to investigate 
the impact of shocks to markets. 

3.4 Granger causality 

Granger causality of a time series {Y t }t=i,...,T with respect to another time series {X t }t=i,...,T 
means that {Yt} provides statistically significant information about {Xt}, in other words, it is 
helpful to predict {X t }. This is often denoted as {Y t } — > {X t } ("{Yt} Granger causes {Xt}"). 
On the other hand, {Y t }t=i,...,T does not Granger cause {X t }t=i,...,T if it does not have any 
explanatory power with respect to future observations Xt+ s , s > 1. 

Granger causality of one time series on another can easily be investigated using the COPAR 
model. As discussed above, the COPAR model directly gives the conditional distribution of 
Xt+i\Xi, ..., Xt, Yi, ..., Yt- If {Yt} does not Granger cause {Xt}, then this distribution is equal 
to that of Xt+i\X\, ...,Xt, which means that all the pairs specifying between-series dependence 



of {Xt} and {Yt} (see Expressions (3.2) and (3.3)) are independent. The model then collapses 



to two independent serial D- vines for {Xt} and {Yt}, respectively. When pair-copulas of these 
two serial D-vines are chosen as for the pairs in the full model, then the two models are nested 
(most copulas, in particular those that will be considered in this paper, include the independence 
copula as special or boundary case) and a standard likelihood-ratio test can be used to investigate 
Granger causality. Let Irv denote the log likelihood of the joint model of {Xt} and {Yt} and 
£dv the cumulative one of the two separate D-vines for {Xt} and {Yt} constructed with the 
same pair-copulas as the joint model, then 

Wrv-Zdv)~4 rv - Pdv , (3.11) 

where Xq denotes a \ 2 distribution with q degrees of freedom and p^y & n d Pdv denote the 
number of parameters of the full (R-vine) and the reduced (two D-vines) model, respectively. 
If copula families are however chosen differently, then a test for non-nested hypotheses such as 



the one by Vuong (1989) can be used. 

To summarize, in order to investigate Granger causality of a time series {^}t=i,....r on 
another time series {Xt}t=i,....T, the COPAR model can be used in conjunction with a likelihood- 
ratio test. If Granger causality of {Xt} on {Yt} is to be investigated, the roles of the two time 
series have to be interchanged. 

4 Applications 

In this section we discuss three relevant applications. First, we analyze four monthly macro- 
economic indicators pairwisely, namely inflation and interest rates as well as stock returns and 
industrial production. Since this is the classical area of application of VAR models it is partic- 
ularly interesting to see what COPAR models can add here. Second, COPAR models are used 
to model daily electricity load demands in four Australian states. Due to geography there is a 
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strong interdependence among states which needs to be captured. Finally, monthly Fama bond 
portfolio returns with medium duration are analyzed. 

In each of the three data sets we compare the COPAR model to relevant benchmark mod- 
els in terms of out-of-sample predictive ability: First, VAR(A;) models are fitted using the 
-R-package vars (iPfaff 2008). Second, we also fit standard copula models with AR(fe) and 



AR(fc)-GARCH(l,l) margins, where the distribution of the innovations is chosen as in the 
corresponding COPAR model and the copula used is selected according to the AIC from a 
range of different copulas capturing all types of dependence (tail-independent Gaussian (N) and 
Frank (F), symmetric-tail-dependent Student-t, lower-tail-dependent Clayton (C), upper-tail- 
dependent Gumbel (G) and Joe (J), and their survival versions (SC, SG, SJ)). Note that the 
copula model with AR(£;)-GARCH(1,1) is a very tough competitor model, since it allows for 
time- varying variances, which the COPAR and the VAR models do not. Since the use of copula- 
GARCH models is a common tool for the analysis of financial time series, it is also included in 
the analysis. 

For COPAR models we distinguish 

• unconditional prediction of Xt+i given Xi, ...,Xt, Y\, ..., Yr, 

• joint prediction of Xt+i and lr+i given X\, ..., Xt, Y\, ..., Yr, and 

• conditional prediction of lr+i given Xi, ..., Xt+i, Y\, ..., Yt- 

To speed up computations IFM-type sequential estimation as described in Section |3.2| is used 
for out-of-sample prediction, since the model is re-estimated for each additional prediction. In 
all three applications model parameters estimated in this way were very close to full maximum 



likelihood estimates. Copulas are selected as described in Section 3.2 and from the same list as 
above. 

Predictions are evaluated in terms of two different loss functions. On the one hand, we 
consider the classical root mean squared error 



RMSE(a;;x) 



N 



T* 

2 
i) > 



1 
7""* / v V. 



rr-\ ■ ' \*i •& 

4=1 



where T* denotes the number of out-of-sample predictions and x = (x±, ...,xt*)' are point 
forecasts of x = (x\, ...,xt*)'- On the other hand, in order to check the coverage of empirical 



prediction intervals, mean interval scores by Gneiting and Raftery (2007) are taken into account: 

T* r 

1 . 2 2 

MlS a (l,u;x) = — Y, («i - U) +-{U ~ Xi)l {x . <t . } + -{*i- Ui)l { u i<Xi} 



t=x 



a i*^'*' a 



where I = (li, ..., It*)' and u = (ux, ...,ut*)', and [k,Ui], i = 1, ...,T*, are 100(1 — a) % prediction 
intervals. For the copula-based models these are determined as (a/2) and (1 — a/2) sample 
quantiles; for the VAR model closed form expressions using normal quantiles are used. 
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Figure 5: Macro-economic indicators: Time series. 

4.1 Macro-economic indicators 

The first application investigates the interdependencies of four monthly US macro-economic 
indicators: continuously compounded real returns on the S&P500 index (MKT), real interest 
rates of 30-day US Treasury Bills (RRF), continuously compounded growth rates of the US CPI 
(INF) and continuously compounded growth rates of US industrial production (IPG). The data is 
available at http://www.economagic.com/ and covers 493 monthly observations from December 
1959 to December 2000, which we split into a training set of 393 and a testing set of 100 



observations. Analyses of comparable data can be found in Lee (1992) and in Zivot and Wang 



(2006, Chapter 11) 



The time series are shown in Figure [5j From the very different behavior of the series one 
may expect only little dependence between them. A look at rank correlations (not corrected 
for serial dependence) confirms this, since most empirical values are below 0.1 in absolute terms 
and only between inflation rates and stock returns as well as inflation rates and interest rates 
there is a somewhat stronger negative dependence of —0.21 and —0.26, respectively. Therefore, 
it will be particularly interesting to analyze Granger causality for this data. 

For the marginal time series we choose skew-normal distributions to account for skewness in 
the data, a common feature of many economic time series. In particular, we use the parametriza- 



tion by Azzalini (1985) , where a shape parameter controls the skewness. If the shape parameter 
is 0, there is no skewness and the distribution reduces to the normal one. 

In the following, we investigate the causal relationships of all twelve (ordered) pairs of these 
four variables to see how the indicators influence each other. 

We select the lag length as described in Section |3.2| Three of the twelve diagnostic plots are 
shown in Figure [6j To facilitate estimation, all pair-copulas were chosen as Gaussian; Gaussian 
pair-copulas are chosen as proxy, since they can model positive as well as negative symmetric 
dependence and only use one parameter. The plots indicate an appropriate autoregressive order 
of k = 2. 

We then fitted a COPAR(2) model for the training data as described in the introduction to 
this section. Selected copula types and estimated parameters are shown in Table [7] in Appendix 
[BJ Although many tail-symmetric copulas (Gaussian, Student-t, Frank) were chosen, this de- 
pendence is mostly non-Gaussian and there is also evidence for some asymmetric dependence, 



e.g. the serial dependence of inflation rates. AIC, BIC and HQC (3.10) confirmed this: Us 



ing non-Gaussian copulas clearly improved the in-sample fit over an only Gaussian COPAR(2) 
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Figure 6: Macro-economic indicators: Information criteria (3.10) for the fit of three variable 



pairs at different lag lengths, where all pair-copulas are chosen as Gaussian. 
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Figure 7: Macro-economic indicators: Granger causality of variables according to the COPAR(2) 
model (left panel) and the VAR(2) model (right panel) at the 5% level. 

model. In addition, there is also clear indication of skewness in the data as determined by the 
estimated marginal parameters. 

Next we compare unconditional and joint prediction using the COPAR(2) model to pre- 
dictions obtained from VAR(2), copula-AR(2) and copula-AR(2)-GARCH(l,l) models on the 
testing set of 100 out-of-sample observations. In addition, we investigate whether predictions 
can be improved through conditional prediction, when the future realization of one variable is 
already known. 

Tables [T] and [2] show root mean squared errors and mean interval scores for all models and all 
twelve variable pairs based on lOO'OOO samples. While the COPAR(2) model almost exclusively 
outperforms the VAR(2) model, the inclusion of GARCH(1,1) effects in the standard copula- 
AR(2) model still provides better predictions. Nevertheless, the standard copula-AR(2) model 
is inferior to the COPAR(2) model which indicates that serial dependence is not linear unlike 
modeled by the AR-margins. Rather surprisingly, conditional prediction is not useful for this 
data. This however makes sense in light of the small dependence between series described above. 

Finally, we investigate Granger causality between the time series using the likelihood ratio 



test (3.11). The resulting causalities as identified by COPAR(2) and VAR(2) models are in- 



dicated by arrows in Figure [7| The two models mainly identify the same interdependencies, 
differences are likely due to the non-linear modeling of the COPAR model. In particular, infla- 
tion rates are determined to significantly influence all other variables, while industrial production 
growth only causes interest rates. Nonetheless, although there is Granger causality among most 
variable pairs, conditional predictions did not improve the forecasting accuracy. 
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MKT|RRF 


MKT|INF 


MKT|IPG 


RRF|MKT 


RRF|INF 


RRF|IPG 


uncond. COPAR(2) 


0.02908 


0.02869 


0.02920 


0.00043 


0.00041 


0.00042 


joint COPAR(2) 


0.02924 


0.02861 


0.02921 


0.00044 


0.00040 


0.00042 


VAR(2) 


0.02900 


0.02867 


0.02923 


0.00042 


0.00042 


0.00043 


copula-AR(2) 


0.02943 


0.02948 


0.02945 


0.00046 


0.00045 


0.00046 


copula-AR(2)-GARCH(l,l) 


0.02897 


0.02893 


0.02895 


0.00041 


0.00041 


0.00041 


cond. COPAR(2) 


0.02924 


0.02912 


0.02928 


0.00043 


0.00041 


0.00043 





INF|MKT 


INF|RRF 


INF|IPG 


IPG|MKT 


IPG|RRF 


IPG|INF 


uncond. COPAR(2) 


0.00144 


0.00141 


0.00144 


0.00437 


0.00441 


0.00443 


joint COPAR(2) 


0.00144 


0.00143 


0.00143 


0.00433 


0.00437 


0.00443 


VAR(2) 


0.00148 


0.00144 


0.00147 


0.00446 


0.00445 


0.00445 


copula-AR(2) 


0.00156 


0.00156 


0.00156 


0.00498 


0.00498 


0.00498 


copula- AR(2)-GARCH(1,1) 


0.00145 


0.00146 


0.00146 


0.00437 


0.00438 


0.00437 


cond. COPAR(2) 


0.00152 


0.00141 


0.00143 


0.00437 


0.00434 


0.00444 



Table 1: Macro-economic indicators: Root mean squared errors of unconditional, joint and 
conditional predictions from COPAR(2) models as well as of predictions from VAR(2) and 
copula- AR(2)-(GARCH(l,l)-)models for all twelve variable pairs and 100 out-of-sample values. 
Best performing methods (other than the conditional) are indicated in bold. 





MKT|RRF 


MKT|INF 


MKT|IPG 


RRF|MKT 


RRF|INF 


RRF|IPG 


uncond. COPAR(2) 


0.1688 


0.1639 


0.1693 


0.0024 


0.0023 


0.0023 


joint COPAR(2) 


0.1707 


0.1649 


0.1705 


0.0026 


0.0027 


0.0028 


VAR(2) 


0.1649 


0.1606 


0.1646 


0.0029 


0.0029 


0.0029 


copula-AR(2) 


0.1703 


0.1712 


0.1706 


0.0030 


0.0030 


0.0030 


copula- AR(2)-GARCH(1,1) 


0.1611 


0.1611 


0.1622 


0.0021 


0.0021 


0.0021 


cond. COPAR(2) 


0.1716 


0.1626 


0.1705 


0.0026 


0.0026 


0.0026 





INF|MKT 


INF|RRF 


INF|IPG 


IPG|MKT 


IPG|RRF 


IPG|INF 


uncond. COPAR(2) 


0.0084 


0.0083 


0.0084 


0.0312 


0.0316 


0.0308 


joint COPAR(2) 


0.0084 


0.0082 


0.0084 


0.0318 


0.0316 


0.0310 


VAR(2) 


0.0094 


0.0092 


0.0093 


0.0317 


0.0320 


0.0311 


copula- AR( 2) 


0.0094 


0.0093 


0.0094 


0.0325 


0.0325 


0.0325 


copula- AR(2)-GARCH(1,1) 


0.0089 


0.0089 


0.0089 


0.0284 


0.0285 


0.0286 


cond. COPAR(2) 


0.0086 


0.0081 


0.0084 


0.0316 


0.0310 


0.0310 



Table 2: Macro-economic indicators: Mean interval scores of unconditional, joint and conditional 
95% prediction intervals from COPAR(2) models as well as of 95% prediction intervals from 
VAR(2) and copula- AR(2)-(GARCH(l,l)-)models for all twelve variable pairs and 100 out-of- 
sample values. Best performing methods (other than the conditional) are indicated in bold. 
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Figure 8: Electricity load demands: Time series. 





QLD 


NSW 


VIC 


SA 


QLD 


1.00 


0.40 


0.13 


0.08 


NSW 


0.40 


1.00 


0.49 


0.39 


VIC 


0.13 


0.49 


1.00 


0.68 


SA 


0.08 


0.39 


0.68 


1.00 



Table 3: Electricity load demands: Unconditional empirical rank correlations between time 
series (not corrected for serial dependence). 



4.2 Australian electricity load demands 

In this application we analyze average daily electricity load demands in Gigawatt for the Aus- 
tralian states Queensland (QLD), New South Wales (NSW), Victoria (VIC) and South Australia 
(SA). The observations have been calculated by averaging the half- hourly observed data for one 
day which are available at http://www.aemo.com.au/. The observed time period is from May 
16, 2005 to June 30, 2008 with 1135 daily observations, split into a training set of 635 and a 
testing set of 500 observations. Figure [8] shows the time series after preprocessing the data to 
remove trend and weekly and annual seasonalities. 

To allow for possible skewness in the data we again choose the skew-normal distribution by 



Azzalini (1985) for the margins. In contrast to the previous application there however is clear 
between-series dependence. Unconditional empirical between-series dependence (not corrected 
for serial dependence in the margins) is shown in Table [31 Since the Australian infrastructure is 
concentrated near the coast, dependence is observed in particular for the pairs QLD-NSW, NSW-VIC 
and VIC-SA, which mirror the geographical location of the states on the coastline. While focusing 
on these pairs, we investigate the causal relationships of all twelve ordered pairs of these four 
states. 

The lag length of COPAR(A;) models for daily electricity load demands is selected as k = 2 



according to different information criteria as described in Section 3.2 The comparison of CO- 
PAR(2) models with Gaussian copulas and with arbitrary copulas then shows that dependence 
is mainly tail-symmetric but mostly non-Gaussian. In particular, tail-dependent Student-t cop- 
ulas play a central role in describing the between-series dependence. Such tail dependence is 
ignored by Gaussian copulas and therefore again stresses the need for non-Gaussian modeling. 
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QLD|NSW 


QLD|VIC 


QLD|SA 


NSW|QLD 


NSW|VIC 


NSW|SA 


uncond. COPAR(2) 
joint COPAR(2) 
VAR(2) 


0.653 

0.663 
0.674 


0.648 

0.670 
0.670 


0.645 

0.657 
0.668 


0.692 

0.706 
0.697 


0.674 

0.674 
0.690 


0.679 

0.689 
0.690 


copula- AR(2) 

copula- AR(2)-GARCH(1,1) 


0.673 
0.679 


0.674 
0.679 


0.673 
0.679 


0.697 
0.699 


0.697 
0.699 


0.697 
0.698 


cond. COPAR(2) 


0.606 


0.624 


0.627 


0.654 


0.603 


0.639 





VIC|QLD 


VIC|NSW 


VIC|SA 


SA|QLD 


SA|NSW 


SA|VIC 


uncond. COPAR(2) 


0.676 


0.674 


0.666 


0.691 


0.686 


0.691 


joint COPAR(2) 


0.672 


0.679 


0.684 


0.690 


0.686 


0.679 


VAR(2) 


0.693 


0.692 


0.685 


0.687 


0.684 


0.688 


copula- AR(2) 


0.690 


0.690 


0.689 


0.686 


0.687 


0.685 


copula-AR(2)-GARCH(l,l) 


0.692 


0.692 


0.692 


0.687 


0.687 


0.687 


cond. COPAR(2) 


0.638 


0.598 


0.555 


0.666 


0.644 


0.560 



Table 4: Electricity load demands: Root mean squared errors of unconditional, joint and con- 
ditional predictions from COPAR(2) models as well as of predictions from VAR(2) and copula- 
AR(2)-(GARCH(l,l)-)models for all twelve variable pairs and 500 out-of-sample values. Best 
performing methods (other than the conditional) are indicated in bold. 



Selected copula types and estimated parameters are shown in Table [8] in Appendix [Bj 

The out-of-sample predictive ability of the COPAR(2) model is evaluated on the testing set 
of 500 observations. Evaluation criteria are computed based on lOO'OOO samples. Tables [4] and 
p| show that the COPAR(2) model performs very strongly and clearly outperforms the VAR(2) 
model and in most cases even the copula- AR(2)-GARCH(1,1) model. 

Due to the clear geographical relationships between the states, conditional prediction clearly 
improves the forecasting accuracy. In particular when the demand of the neighboring state is 
known, there is a significant improvement. This result indicates that the COPAR model may 
be useful for scenario analysis, e.g. to examine the effect of shocks like an extreme demand to 
the electricity market. 

4.3 Fama bond portfolio returns 

After analyzing economic time series data in first two applications, we now investigate the inter- 
dependencies among financial variables, namely monthly Fama bond portfolio^ returns obtained 
from the CRSP Center for Research in Security Prices for maturities 12 and 24 months (variables 
12M and 24M, respectively). The time series are observed from January 1952 to December 2003 
with 684 monthly observations (training set of 584, testing set of 100 observations) and shown 
in Figure M 



Artificial zero-coupon bonds constructed after first extracting the term structure from a subset of US bonds. 
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QLD|NSW 


QLD|VIC 


QLD|SA 


NSW|QLD 


NSW|VIC 


NSW|SA 


uncond. COPAR(2) 


3.567 


3.551 


3.540 


3.771 


3.838 


3.836 


joint COPAR(2) 


3.701 


3.705 


3.612 


3.818 


3.790 


3.887 


VAR(2) 


3.901 


3.855 


3.839 


3.788 


3.987 


3.968 


copula- AR(2) 


3.881 


3.883 


3.865 


3.789 


3.797 


3.792 


copula- AR(2)-GARCH(1,1) 


3.697 


3.706 


3.705 


3.843 


3.849 


3.845 


cond. COPAR(2) 


3.126 


3.291 


3.258 


3.292 


3.039 


3.280 





VIC|QLD 


VIC|NSW 


VIC|SA 


SA|QLD 


SA|NSW 


SA|VIC 


uncond. COPAR(2) 


3.766 


3.781 


3.834 


3.262 


3.237 


3.288 


joint COPAR(2) 


3.770 


3.809 


3.881 


3.299 


3.287 


3.235 


VAR(2) 


3.806 


3.782 


3.847 


3.275 


3.236 


3.280 


copula- AR(2) 


3.817 


3.798 


3.813 


3.248 


3.239 


3.244 


copula-AR(2)-GARCH(l,l) 


3.822 


3.818 


3.822 


3.263 


3.244 


3.252 


cond. COPAR(2) 


3.300 


3.175 


3.074 


3.010 


2.946 


2.658 



Table 5: Electricity load demands: Mean interval scores of unconditional, joint and conditional 
95% prediction intervals from COPAR(2) models as well as of 95% prediction intervals from 
VAR(2) and copula-AR(2)-(GARCH(l,l)-)models for all tewlve variable pairs and 500 out-of- 
sample values. Best performing methods (other than the conditional) are indicated in bold. 



In contrast to the first two applications (low and medium between-series dependence, re- 
spectively), unconditional empirical between-series rank correlation (not corrected for serial de- 
pendence in the margins) between the Fama bond portfolio return series is quite strong (0.84). 

Also in contrast to the first two applications we opt for a more sophisticated marginal 
distribution to appropriately capture the behavior of the financial time series: We choose the 
hyperbolic distribution which is widely used for financial return data and accounts for both fat 



tails and skewness (see McNeil et al. (2005) ). As for the skew- normal distribution, the hyperbolic 



distribution is symmetric if the skewness parameter is 0. An additional shape parameter controls 
how much weight is assigned to the tails and to the center of the distribution. 

Serial dependence in the bond return data is more persistent than in the previous two 
applications. Based on the discussion in Section 3.2 we select a lag length of k = 3 and there 



even is some indication that maybe a higher order may be reasonable. 

The in-sample fit and copula selection shows that there is clear non-Gaussian dependence in 



the data (see Tables^] and 10 in Appendix pi). While Student-t and Frank copulas are central for 
serial and between-series dependence modeling, also tail- asymmetric Joe, Gumbel and Clayton 
copulas have been identified to determine the joint behavior of the time series. 

Prediction evaluation criteria based on the out-of-sample testing set of 100 observations and 
on lOO'OOO samples are shown in Table pi The COPAR(3) model's predictive ability is clearly 
superior to that of the VAR(3) model as well as to a standard copula- AR(3) model with and 
without GARCH(1,1) effects. Also conditional prediction improves the forecasting accuracy due 
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Figure 9: Bond returns: Time series. 





RMSE 


MIS 




12M|24M 


24M|12M 


12M|24M 


24M|12M 


uncond. COPAR(3) 


0.00140 


0.00431 


0.00796 


0.02392 


joint COPAR(3) 


0.00144 


0.00420 


0.00803 


0.02110 


VAR(3) 


0.00165 


0.00457 


0.01069 


0.02595 


copula- AR(3) 


0.00171 


0.00520 


0.01087 


0.02697 


copula-AR(3)-GARCH(l,l) 


0.00372 


0.00472 


0.04986 


0.02344 


cond. COPAR(3) 


0.00081 


0.00230 


0.00411 


0.01207 



Table 6: Bond returns: Root mean squared errors and mean interval scores of unconditional, 
joint and conditional predictions and 95% prediction intervals from COPAR(3) models as well 
as of predictions and 95% prediction intervals from VAR(3) and copula-AR(3)-(GARCH(l,l)- 
)models for both ordered variable pairs and 100 out-of-sample values. Best performing methods 
(other than the conditional) are indicated in bold. 



to the significant between-series dependence. As before this result may be used for scenario 
analyses of the bond market: How do shocks to bonds with one maturity influence bonds with 
the other maturity? 

5 Conclusion 

In this paper we described the novel copula autoregressive model which benefits from the flexi- 
bility of R-vine copulas and allows to model non-linear dependence among multiple time series. 
While three relevant applications to financial and economic time series demonstrate the useful- 
ness of these models, their practical importance will be even more pronounced when stronger 
non- linear and non-symmetric dependencies are present, in time as well as between series. 

The applications however also showed that time-varying variance effects may be needed in 
certain cases and improve over time-constant modeling. Approaches to this issue are subject of 
future research. 
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A Technical supplement 



As noted in Section |3.2| sequential copula selection and likelihood computation is not straight- 
forward in the COPAR model. Figure [4] illustrates the interdependencies among the copulas. 
For likelihood computation, of course the R-vine matrix specification (3.8) could be used in- 
stead. However, given that the number of time points T might be large, evaluation of this 
2T-by-2T-matrix is computationally rather inefficient, since most matrix entries do not contain 
any information due to the assumed autoregressive order (see Example 3.7). 

In the following we therefore present how to sequentially select copulas in a COPAR(2) model 
for data {xt}t=i,...,T and {yt}t=i,...,T- Rather than selecting copulas, preliminarily determined 
copulas can be estimated or likelihoods can be evaluated by simply altering the respective lines. 

(i) Serial dependence of {X t }t=i,...,T'- X s , X t \X s+ i, ..., X t -i, 1 < s < t < T. 

(a) Select copula C* = C Xt _ lXt based on {F x (x t )}t=i,...,T-i and {F(x t )}t=2,...,T- 

(b) Compute F(xt\xt-i), t = 3, ...,T, and F(xt-i\xt), t = 2, ..., T, from C^ using Ex- 
pression \2A\ . 

(c) Select copula C- = C x t -iX t \Xt-i based on F{xt-i\xt), t = 2, ...,T— 1, and F(xt\xt-i), 
i = 3,...,T. 

(d) Compute F(x t -2\x t -i,Xt), t = 3, ...,T, and F(x t \x t -2, x t -i), t = 3,...,T, from Ca- 
using Expression (2.4). 

(ii) Between-series dependence X s ,Y t \X s+ i, ...,X t , 1 < s < t < T. 

(a) Select copula Cq Y = C Xt Y t based on {F x (x t )}t=i,...,T and {F Y (y t )}t=i,...,T- 

(b) Compute F(y t \xt), t = 2, ...,T, from Cq Y using Expression \2A\ . 

(c) Select copula C^ Y = C Xt _ 1 y t \ Xt based on F(xt-i\xt), t = 2, ...,T, and F(yt\xt), t = 
2,...,T. 

(d) Compute F(y t \xt-i,x t ), t = 3, ...,T, from C^ Y using Expression (2.4). 

(e) Select copula C£- Y = C Xt _ 2 Y t \x (t _ 1) . t based on F(xt-2\xt-i,xt), t = 3,...,T, and 
F(y t \x t -i,x t ), t = 3, ...,T. 



(f) Compute F(y t \x t -2, x t -i, x t ), t = 3, ...,T— 1, from C 2 Y using Expression (2.4). 

(hi) Between-series dependence Y s , Xt\X\, ..., Xt-i, Y s +i, ...,Yt-i, 1 < s < t < T, and condi- 
tional serial dependence of {lt}t=i . t' Y a ,Yt\Xi,...,Xt,Y s +i,...,Yt-.i, 1 < s < t < T. 



■iYX 

F(x t \xi,...,x t -i), t = 2,...,T, where 



(a) Select copula C YX = C Yt _ lXt \x 1 . (t _ 1) based on F(y t \xi, ...,xt), t = 1,...,T — 1, and 



F(y t \xi,...,x t ) 



F(yi\x x ) t = \ 

F(y2\xi,x 2 ) t = 2 

F(y t \xt-2,xt-i,x t ) t>3 
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and 



F(x t \xi,...,x t -i, 



F(x 2 \xi] 



t = 2 



F(xt\xt-2,xt-i) t>3 



(b) Compute F(x t \xi, ...,xt-i,yt-i), t = 3, ...,T, and F(y t -i\xi, ...,x t ), t = 2, ...,T, from 
C{ using Expression (2.4). 



(c) Select copula C\ = C Yt _ xYt \x x .. t based on F{y t -\\x\, ...,x t ), t = 2, ..., T, and F(y t \xi, ...,x t ), 
t = 2, ...,T, where 



F(yi|xi,x 2 ) i = 2 

F(y 2 |xi, x 2 ,x 3 ) t = 3- 

F(y t -i\xt-3,...,x t ) t>4 



F(y t _i|xi,...,x t ) : 
(d) Compute F(yj_i|xi,...,x t ,yt), i = 2, ...,T - 1, and F(y t \xi, ...,x t ,yt-i), t = 3,...,T, 



from Cj using Expression (2.4). 

(e) Select copula C\ x = C Yt _ 2 Xt\x 1:(t _ 1} Yt-i based on F(y t -i\xi,...,xt,yt), t = 2,...,T- 
1, and F(xt\xi,...,xt-i,yt-i), t = 3, ...,T, where 



F(y 4 _i|xi,...,x t ,y t ) = < 



^(2/1 \xi,x2,y2) t = 2 

F(y 2 \xi,X2,X3,y 3 ) t = 3, 

F(y t -i\x t -3,xt-2,xt-i,xt,yt) t>4 



and 



F(x t \xi, ...,x t -i,yt-t) = < 



F(x 3 \xi,x 2 ,y2) t = 3 

F(x4\xi,X2,X3,y3) t = A. 

^F(x t \x t -4,Xt-3,Xt-2,Xt-l,yt-l) t > 5 

(f) Compute F(y t -2\xi, ...,xt,yt-i), t = 3, ...,T, from C\ x using Expression (2.4). 



(g) Select copula C\ = Cy t _ 2 y t |x 1:t y t _i based on F(y t _ 2 |xi,...,x t ,2/t_i), i = 3, ..., T, and 
F(y t |xi,...,a;t,yt_i), i = 3, ...,T, where 



F(y t _ 2 |xi,...,x t ,yi_i) 



F(yi|xi,x 2 ,x 3 ,y 2 ) £ = 3 

■^(2/2|a:i 3 2:2,^3, sc 3 , S4,y 3 ) t = 4 

F(y i _ 2 |x t _4,...,x t ,y i _i) £>5 



and 



F(y t |xi,...,xt,y t _i) : 
For autoregressive orders higher than 2 the above procedure can be easily extended. 



F(y3\xi,X2,x 3 ,y 2 ) t = 3 

F(y 4 |xi,x 2 ,X3,X4,y 3 ) t = 4 
F(y t \xt-A,...,x u yt-i) t>5 
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B Additional figures and tables 
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Figure 10: R-vine for X\, ..., X4, Y±, ..., Y4 extended to include X§ and I5. Dashed lines indicate 
additional edges to include X5, dotted those to include I5. 
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F 


t 


t 


t 


RRF 


0.0028 


0.0027 


-1.3316 


1.95 


0.26 
12.18 


-0.09 


-0.06 
30 


-0.02 


0.01 


-0.09 

17.43 


0.93 
4.61 


0.16 
4.62 










049 


0.17 


-0.01 


-0.04 


-0.01 


0.00 


-0.06 


0.77 


0.10 


INF, 


6e-04 


0.005 


2.679 


G 


t 


N 


N 


t 


F 


F 


t 


F 


IPG 


-0.0016 


0.0104 


0.7242 


1.92 


0.26 

12.48 


-0.09 


-0.10 


-0.09 
30.00 


0.14 


-0.06 


0.39 

9.88 


0.93 










0.48 


0.17 


-0.05 


-0.06 


-0.06 


0.02 


-0.01 


0.25 


0.10 


IPG, 


-0.0016 


0.0105 


0.7242 


t 


F 


C 


N 


N 


SG 


SG 


N 


F 


MKT 


0.0353 


0.0488 


-1.8383 


0.39 
9.25 


1.03 


0.04 


-0.09 


-0.10 


1.03 


1.09 


0.25 


-0.49 










0.26 


0.11 


0.02 


-0.06 


-0.06 


0.02 


0.08 


0.16 


-0.05 


IPG, 


-0.0021 


0.0104 


0.7242 


t 


F 


C 


SJ 


SJ 


F 


t 


F 


SG 


RRF 


0.0027 


0.0027 


-1.3316 


0.39 

8.48 


1.01 


0.03 


1.01 


1.03 


-0.07 


-0.08 

8.88 


15.16 


1.17 










0.26 


0.11 


0.02 


0.01 


0.02 


-0.01 


-0.05 


0.76 


0.14 


IPG, 


-0.0021 


0.0104 


0.7242 


t 


F 


N 


N 


N 


t 


F 


G 


t 


INF 


5e-04 


0.005 


2.679 


0.41 
10.31 


0.96 


-0.09 


-0.03 


-0.03 


-0.13 
30.00 


-0.71 


1.91 


0.26 
12.17 










0.27 


0.1 1 


-0.06 


-0.02 


-0.02 


-0.08 


-0.08 


0.48 


0.17 



Table 7: Macro-economic indicators: Estimated parameters and chosen copulas of the CO- 
PAR(2) models, where the first time series in each block corresponds to {Xt} and the second to 
{It}. Italic numbers indicate corresponding Kendall's rs. 
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Margins 












Copulas 












loc. 


scale 


skew 


C* 


cf 


r*XY 


C? Y 


or 


cY x 


cV 


C\ 


cl 


QLD, 


0.5347 


1.1087 


-0.6587 


F 


F 


t 


F 


sc 


SC 


F 


F 


N 


NSW 


-0.8134 


1.2549 


1.0882 


7.05 


-0.55 


0.37 
8.63 


-0.20 


0.05 


0.18 


-0.14 


5.63 


-0.07 










0.56 


-0.06 


0.24 


-0.02 


0.02 


0.08 


-0.02 


0.49 


-0.04 


QLD, 


0.5627 


1.1083 


-0.6831 


F 


F 


SJ 


F 


F 


F 


N 


F 


N 


NSW 


-0.5783 


1.0801 


0.8359 


7.03 


-0.53 


1.14 


0.19 


0.06 


0.01 


0.10 


5.78 


-0.12 










0.56 


-0.06 


0.07 


0.02 


0.01 


0.00 


0.06 


0.50 


-0.08 


QLD, 


0.647 


1.1597 


-0.8191 


F 


F 


SJ 


N 


F 


F 


F 


G 


F 


SA 


-0.9021 


1.4215 


1.9323 


7.05 


-0.54 


1.06 


-0.02 


-0.03 


0.01 


0.62 


1.85 


-0.53 










0.56 


-0.06 


0.04 


-0.01 


0.00 


0.00 


0.07 


0.46 


-0.06 


NSW, 


-0.8824 


1.315 


1.1642 


F 


F 


t 


t 


G 


F 


G 


t 


F 


QLD 


0.7313 


1.2351 


-1.1263 


6.46 


-0.43 


0.39 

9.47 


0.14 
18.29 


1.05 


0.06 


1.04 


0.72 
13.81 


-0.52 










0.54 


-0.05 


0.26 


0.09 


0.05 


0.01 


0.04 


0.51 


-0.06 


NSW, 


-0.6892 


1.2084 


0.8284 


F 


F 


t 


N 


F 


F 


t 


F 


F 


VIC 


-0.6653 


1.113 


1.0071 


6.36 


-0.64 


0.44 
6.99 


-0.16 


-0.03 


2.29 


-0.10 
13.73 


4.70 


-0.58 










0.53 


-0.07 


0.29 


-0.1 


0.00 


0.24 


-0.07 


0.44 


-0.06 


NSW, 


-0.913 


1.3292 


1.2318 


F 


F 


t 


F 


N 


N 


N 


t 


F 


SA 


-0.6843 


1.1507 


1.0412 


6.43 


-0.50 


0.43 
5.57 


-0.51 


0.01 


0.35 


-0.07 


0.56 
8.65 


-0.55 










0.54 


-0.05 


0.28 


-0.06 


0.01 


0.22 


-0.04 


0.38 


-0.06 


VIC, 


-0.6205 


1.0997 


0.9191 


F 


N 


SJ 


F 


N 


F 


C 


t 


F 


QLD 


0.7381 


1.2083 


-1.0317 


5.88 


-0.11 


1.17 


-0.06 


0.04 


0.33 


0.03 


0.76 
9.18 


-0.64 










0.51 


-0.07 


0.09 


-0.01 


0.03 


0.04 


0.02 


0.55 


-0.07 


VIC, 


-0.5819 


1.0952 


0.8255 


F 


N 


t 


t 


F 


N 


F 


F 


F 


NSW 


-0.7386 


1.2001 


1.0095 


6.11 


-0.12 


0.43 
8.66 


0.31 
8.33 


0.47 


-0.07 


0.83 


5.17 


0.01 










0.52 


-0.07 


0.28 


0.20 


0.05 


-0.05 


0.09 


0.47 


0.00 


VIC, 


-0.6299 


1.1162 


0.9363 


F 


N 


t 


N 


F 


t 


N 


G 


t 


SA 


-0.7868 


1.2349 


1.4554 


5.87 


-0.11 


0.67 
24.46 


-0.09 


0.26 


0.23 
14.38 


-0.12 


1.48 


-0.03 

14.37 










0.51 


-0.07 


0.47 


-0.06 


0.03 


0.15 


-0.08 


0.32 


-0.02 


SA, 


-0.526 


1.0738 


0.8463 


t 


N 


SJ 


F 


F 


F 


SC 


t 


F 


QLD 


0.7778 


1.2763 


-1.1955 


0.63 
10.76 


-0.07 


1.07 


0.01 


0.37 


0.00 


0.09 


0.76 
10.11 


-0.56 










0.43 


-0.05 


0.04 


0.00 


0.04 


0.00 


0.04 


0.55 


-0.06 


SA, 


-0.6242 


1.1364 


1.0308 


t 


N 


t 


N 


F 


N 


N 


F 


C 


NSW 


-0.758 


1.2212 


1.079 


0.63 
10.76 


-0.07 


0.30 
5.30 


0.31 


0.72 


-0.04 


0.08 


5.28 


0.01 










0.44 


-0.05 


0.19 


0.20 


0.08 


-0.03 


0.05 


0.47 


0.01 


SA, 


-0.5563 


1.0947 


0.8318 


t 


N 


t 


t 


F 


N 


SC 


F 


N 


VIC 


-0.6184 


1.0878 


0.9467 


0.64 
10.53 


-0.07 


0.65 

24.75 


0.26 
8.63 


0.21 


-0.04 


0.05 


4.47 


-0.04 










0.44 


-0.05 


0.45 


0.17 


0.02 


-0.02 


0.02 


0.42 


-0.02 



Table 8: Electricity load demands: Estimated parameters and 
models, where the first time series in each block corresponds 
Italic numbers indicate corresponding Kendall's rs. 
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chosen copulas of the COPAR(2) 
to {Xt} and the second to {!*}. 





shape 


loc. 


scale 


skew 


M12, 


0.6872 


0.0029 


0.0032 


0.0021 


M24 


0.0368 


0.0027 


0.0063 


0.0029 


M24, 


0.0332 


0.0024 


0.0060 


0.0027 


M12 


0.6864 


0.0025 


0.0027 


0.0023 



Table 9: Bond returns: Estimated marginal parameters of the COPAR(3) models. 





c x 


c? 


c? 


r<XY 


C XY 


c? Y 


cf Y 


c YX 


C YX 


cV 


C Y 


C Y 


c Y 


M12, 


F 


t 


F 


G 


t 


F 


t 


t 


t 


F 


t 


F 


SG 


M24 


5.32 


0.32 

5.55 


2.07 


2.90 


-0.27 
13.96 


-1.24 


-0.08 
22.13 


-0.36 
10.51 


-0.27 
18.47 


-0.93 


0.51 
12.96 


1.76 


1.22 




048 


0.21 


0.22 


0.66 


-0.17 


-0.14 


-0.05 


-0.24 


-0.18 


-0.10 


0.34 


0.19 


0.18 


M24, 


J 


t 


F 


G 


SC 


t 


t 


SC 


SC 


SC 


t 


t 


t 


M12 


1.34 


0.09 
8.09 


1.32 


2.59 


0.15 


0.09 
18.07 


0.07 

17.87 


0.35 


0.12 


0.07 


0.76 
2.99 


0.43 
16.98 


0.34 
3.37 




0.16 


0.06 


0.14 


0.61 


0.07 


0.06 


0.04 


0.15 


0.05 


0.03 


0.55 


0.28 


0.22 



Table 10: Bond returns: Estimated parameters and chosen copulas of the COPAR(3) models, 
where the first time series in each block corresponds to {Xt} and the second to {Yt}. Italic 
numbers indicate corresponding Kendall's rs. 
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